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THE FAR-INFRARED, UV AND MOLECULAR GAS RELATION IN GALAXIES UP TO Z=2.5. 



(N 

o 

(N 

> 
O 

in 

(N 

O 

u 

Of 

6 



> 
o 

in 



(N 



R. Nordon 1 , D. Lutz 2 , A. Saintonge 2 , S. Berta 2 , S. Wuyts 2 , N. M. Forster Schreiber 2 , R. Genzel 2 , B. Magnelli 2 , A. Poglitsch 2 , 

P. Popesso 2 , D. Rosario 2 , E. Sturm 2 , L. J. Tacconi 2 
Draft November 27, 2012 

ABSTRACT 

We use the infrared excess (IRX) FIR/UV luminosity ratio to study the relation between the effective UV 
attenuation (Airx) and the UV spectral slope (/?) in a sample of 450 1 < z < 2.5 galaxies. The FIR data is from 
very deep Herschel observations in the GOODS fields that allow us to detect galaxies with SFRs typical of 
galaxies with log(M»)>9.3. Thus, we are able to study galaxies on and even below the main SFR-stellar mass 
relation (main sequence). We find that main sequence galaxies form a tight sequence in the IRX-/? plane, which 
has a flatter slope than commonly used relations. This slope favors a SMC-like UV extinction curve, though 
the interpretation is model dependent. The scatter in the Airx-/? plane, correlates with the position of the 
galaxies in the SFR-M, plane. Using a smaller sample of galaxies with CO gas masses, we study the relation 
between the UV attenuation and the molecular gas content. We find a very tight relation between the scatter in 
the IRX-/? plane and the specific attenuation Sa> a quantity that represents the attenuation contributed by the 
molecular gas mass per young star. Sa is sensitive to both the geometrical arrangement of stars and dust, and 
to the compactness of the star forming regions. We use this empirical relation to derive a method for estimating 
molecular gas masses using only widely available integrated rest-frame UV and FIR photometry. The method 
produces gas masses with an accuracy between 0.12-0.16 dex in samples of normal galaxies between z~0 and 
z~1.5. Major mergers and sub-millimeter galaxies follow a different Sa relation. 

Subject headings: Galaxies: evolution - Galaxies: starburst - Galaxies: fundamental parameters - Cosmology: 
observations - Infrared: galaxies - Ultraviolet: galaxies 



1. INTRODUCTION 

The ultra-violet (UV) part in the spectrum of star forming 
galaxies (SFG), between the Lyman edge and ~3000 A, is 
dominated by emission from young stars and provides strong 
constraints on the ongoing star formation rate (SFR). The 
availability of many SFR tracers diminishes at moderately- 
high redshifts (z > 1), where nebular lines are red shifted into 
the near-infrared and spectroscopy becomes increasingly dif- 
ficult, while the availability of deep mid and far infrared (MIR 
and FIR) photometry is limited to a few deep extragalactic 
fields. In contrast, at these redshifts the UV band shifts into 
the optical range, which makes it possible to observe from the 
ground to very high depths, often making it the most accessi- 
ble SFR indicator for large samples. 

UV radiation is susceptible to attenuation by dust in the host 
galaxy. In many examples, only a few percent of the UV pho- 
tons escape the host galaxy. It is therefore vital to correct 
UV luminosities for dust attenuation. The estimation of the 
attenuation normally relies on measuring the change in some 
photometric colors, often at wavelengths other than the ones 
to be corrected. This estimation also requires assumptions on 
the extinction curve and the geometrical arrangement of the 
UV sources and dust. 

In particular, the relation between the UV attenuation and 
the UV spec tral slope /? ( f» <x AP) has been empirica lly cali- 
brated (e.g., iMeurer et alll999tlCdzetti et al.ll20"00l), mostly 
on loc al galaxies, and used in z > 1 studies ("e.g., lDaddietal.1 
2004). The calibrations used the ratio of the FIR luminosity 
(representing the UV radiation reprocessed by dust) to the UV 
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luminosity (escaped light) as an indicator of the UV attenua- 
tion. The FIR/UV ratio is often referred to as the infrared 
excess or IRX and we label the UV attenuation derived from 

the FIR/UV ratio as A IRX - 

Since the works oflMeurer et alj d 19991 M99 hereafter) and 
ICalzetti et all (120001 COO hereafter) on the At RX -/? relation 
with the International Ultraviolet Explorer (IUE), Infrared 
Space Observatory (ISO) and Infrared Astronomical Satellite 
(IRAS) data, a lar ge body of work has been written on t he sub - 
ject, for example : [Chariot & Falil (120001), iKong et al l (2004), 
Buat et al l (120051) . Seibert et alj (120051). Howell et all (hoi 
Hao et alJ <l20TTT). Iwild et all (1201 11) . (Overzier et all 
Shimiz uet alJ (|201 1), and references therein 



2011b . 
The need for 

far infrared luminosities in order to determine Airx signifi- 
cantly limited earlier studies to low redshifts. At redshifts of 
z > 1 such studies were limited to very FIR-bright objects, 
or to use alternative methods for determining the UV atten- 
uation. The situation improved with the advent of Spitzer 
(used by some of the ab ove references) and recently Her- 
schel (Pilb ratt et alJ l2010) infrared space telescopes. The lat- 
ter in particular performed deep extragalactic surveys that de- 
tect normal star forming galaxies (SFGs) at redshifts of 2 and 
more with a wavelength coverage from 70 to 500 /vm in six 
bands that cover the peak of the FIR emission in SFG. Some 
of the recent Hersc hel results regarding t he UV attenuatio n 
and Airx-/? include: iNordon et alJ (120101) iBuat et all (120101 
Wiiesingh e et all d201 ll). iBurgarella et alJ (120111) , IBuat etal 
d2011l) .lReddv et al l (120121) . 

Most of the above mentioned studies report significant scat- 
ter in the Airx-/? relation, with different galaxy populations 
tending to lie on different relations. This may indicate differ- 
ent dust properties, different stellar populations that dominate 
the UV, or different star and dust geometry in these popula- 
tions. The latter is of great importance as the arrangement of 
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stars and dust is not yet clear and may have a strong effect 
on AiRx-y6. Some authors favor a dust screen in the fore- 
ground with the UV emitting stars behind it. Others find 
better agreement with a mix of dust and stars in the same 
volume, or com binations of both possibilities. For example, 
iBuat & Xul (119961) adopt a model where the stars and dust are 
mixed up to some scale height in the galaxy. Above this layer, 
additional dust acts a s a foreground screen. The model of 
ICharlot & Falll (2000) incorporates an age-dependent obscu- 
ration where the younger stars are more obscured than older 
stars that already left, or cleared away the molecular clouds in 
which they formed. 

A relatively tight relation has been found between the stel- 
lar mass of galaxies and their SFR. This relation evolves 
with redshift and has been shown to be in place already at 
z=2 jBrinchmann et al. |[ 2004t iNoeske et al.l l2007: Elba z et all 
120071 iDaddi et alJl2007l) . The relation is often referred to as 
the main-sequence of star forming galaxies. In this paper 
we will refer to it as the 'main sequence' (MS) for brevity. 
Main sequence galaxies at a given redshift show many simi- 
larities to each other, in particular in properties that relate to 
star formation: e.g., they tend to show similar mid- to far- 
infrared SEDs (Elba z et all 1201 U iNordon et all 120 121) SFR 
surface dens ities (Elbaz et al. 2011; Ruiopakarn et al.l 1201 It 



TABLE 1 

PACS PHOTOMETRY DEPTHS IN THE FIELDS USED IN THIS STUDY. 

3<r flux limits [mJy] 
Field 70 pm 100 pm 160 pm 



Wuvts et all 1201 ll), gas depletion time scales ( Daddi et all 
2010t iTacconi et alj|2010t ISaintonge et al.ll201 lbl) . light pro- 



file dWuvts et al.ll201 1 ), etc. These properties tend to change 
for galaxies with an enhanced SFR over the main relation (i.e. 
'above' the main sequence). It has been known that galaxies 
with different 'birthrate' parameters b = S FR/ < S FR > (ra- 
tio of current to past SFR, where < SFR > is the mass over 
time since the initial formation of the galaxy) te nd to lie on 
different A IRX -j0 relations (e.g., iKong et ail20 04). At a given 
redshift, this parameter is closely related to the specific SFR 
(SSFR-S FR/M,). Since main sequence galaxies at a given 
redshift tend to have a narrow range of SSFR, one may ex- 
pect an AiRx-y6 relation for main sequence galaxies. Such a 
relation is expected to have some scatter to it, in particular at 
high redshifts (z > 3) where time scales are short and the star 
formation history c an have a significant effect on the inferred 
physical quantities (Schaerer et al. 2012|). 

Stars form inside molecular clouds and a relatively tight re- 
lation has been found between the molecular gas surface den- 
sity a nd the SFR surface density (Schmidt 1959]; iKennicutH 
1989), often referred to as the Kennicutt-Schmidt (KS) rela- 
tion. Resolved giant molecular clouds (GMCs) in the Milky- 
Way and nearby galaxies show high optical depths and so 
the gas and dust in which the stars form also attenuate their 
UV emission. One therefore may expect a relation between 
the SFR, the molecular gas and the attenuation. However, 
there has been little evidence for a direct relation between 
the molecular gas and the UV attenuation, except for a gen- 
eral broad trend of i ncreasing attenuation with the molecula r 
gas surface density (IBuat & Xulll996t iBoissier et al.l I20O7T) . 
I Wuvts et aTl (120111) used a simple model to predict the FIR 
to UV luminosity ratio from the SFR surface density. The 
model used the KS relation to derive the molecular gas sur- 
face density from the SFR. Then, with some assumptions on 
the geometry, they derived the column density for attenuation. 
Their model was fairly successful at low redshifts, but devel- 
oped increasing deviations as redshift increased. The nature 
of such an attenuation-gas relation should hold some clues 
as to the geometrical arrangement of the star and attenuating 
material, yet it is still unclear. 
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In this study we use FIR photometry from t he Photodetec- 
tor A rray Camera and Spectrometer (PACS, iPoglitsch et alJ 
120101) on board Herschel, combined with UBVIz photometry 
to study the Airx-jS in 1 < z < 2.5 galaxies. When combined 
with the broad and deep optical coverage of these fields, the 
PACS observations at 70, 100 and 160 um, reaching depths 
and resolutions never available before, allow us to build a siz- 
able sample to study the Amx-f3 relation and its scatter. Thus, 
we are able to perform this study at the redshifts where rest- 
UV observations become highly useful and often the most ac- 
cessible SFR indicator. In addition, we study the relation be- 
tween molecular gas content from CO measurements and the 
UV attenuation and spectral slope in 1 < z < 1.5 normal star 
forming galaxies (SFG). We also test these relations on low 
redshift galaxies. 

The paper is structured as follows: In section|2]we describe 
the data and samples used in this study. In section[3]we study 
the Airx-j6 diagram and the structure of the scatter in it. In 
section|4]we find a relation between Airx, /3 and the molecu- 
lar gas mass. Then, in section [5] we derive a method to esti- 
mate molecular gas masses based on FIR and UV luminosi- 
ties, which we test on a low redshift sample. In section [6] we 
discuss the effect of various stars and dust configurations on 
the Airx-jS relation and the reason we can determine gas mass 
from FIR and UV photometry. Our conclusions are summa- 
rized in section [7] 

Thro ughout this paper we assume a Chabrier 
(2003) IMF, and a cosmology with (Q„„Q A ,i/ ) = 
(0.3, 0.7, 70^8-! Mpc- 1 ). 

2. DATA AND SAMPLES 

Our far-infrared data is based on Herschel-PACS ob- 
servations in the GOODS fields. We combine the data 
from two Herschel large programs: the PACS E volutionary 
Probe (PEPl|) project, a guaranteed time program dLutz et alJ 
1201 1 | ) and the G OODS-Herschel open time key program 
(Elba z~eTaT1l201 lb. The c ombin ed reduction is described in 
details in iMagnelli et al.l d2012l) and is similar i n principle 
to the reduction of the PEP data as described in dLutz et al] 
120111) . The limiting fluxes are listed in Table [Q The PACS 
160, 100 and 70 //m fluxes were extracted using sources from 
a Spitzer MI PS 24 //m catalog as p riors, foll owing the method 
described in Magnel lTetail ( 120091) (see also lLutz et"aTl d2011l) 
for more details). For our sample selection we require a red- 
shift (either photometric or spectroscopic), a stellar mass, a 
PACS detection, and sufficient optical bands to derive rest 
frame 1600 A and 2800 A magnitudes. Below we describe 
the constraints each requirement imposes on our sample se- 
lection. 

Pho tometric redshifts are the same as used in lNordon et alJ 
d2012l) and thanks to the large number of bands and the high 
limiting magnitudes, their quality is very good. When com- 
pared to the spectroscopic redshifts: {z p i w t - z spe c)/(l + z, pec ) 
= 0.033 and 0.038 in GOODS-N and GOODS-S respectively. 

2 http: //www . mpe . mpg . de/ir/Research/PEP/ 
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63% and 52% of the selected galaxies in GOODS-S and 
GOODS-N respectively have spectroscopic red shifts. 

The stellar masses are the same as used in iNordon et alj 
(12011. They were calculated using the method described 
in iFontana et al.l (120041) . with a djustments as d escrib ed in 
Santin i et al.1 d2009l) . see also IFontana et al.l (120061) and 
Grazian et al.l (|2006). To ensure good mass estimates, we re- 
quire a detection in all of the 3.6-5.8 fim IRAC bands that 
cover the SED rest frame stellar bump (rest frame H-band) 
in z > 1 galaxies. We place a conservative lower limit on 
the selected galaxies masses of log(M»/M ) > 9.3 in order 
to avoid the galaxies with the lowest derived masses that are 
more likely to be critical errors in the mass fit. This limit en- 
sures detections in most of the bands between U and IRAC 
5.8 yum and thus reliable mass estimates. 

Lre (8-1000 fjm) is calculated by fitting IChary & Elbazl 
(120011 CE01 hereafter) templates to the available PACS 
fluxes. Both template shape and scale are allowed to vary 
and the solution with minimal x 2 value is selected. To derive 
errors on Lir we randomize the observed fluxes according to 
their errors while assuming a normal distribution and repeat 
the template fitting to derive a new luminosity. This step is 
repeated a large number of times and the standard deviation 
(STDEV) of the measured luminosities is adopted as the er- 
ror on Lir. If only one PACS band is available (16% of the 
galaxies, virtually all with only 160 /urn), we fit a CE01 tem- 
plate while preserving the original scaling of the templates. 
This introduces only a small additional e rror and no s ignifi- 
cant bias at redshift s 1.0 < z < 2.5 (se e lElbaz et all (120101) 
and the appendix in lNordon et al.l (120121) 1. 

We estimate rest frame 1600 A and 2800 A magnitudes 
(Mi6O0 an d M2&00 respectively) by interpolating between the 
available bands. We use UBVIz photometry from various in- 
struments. In GOODS-N a photometric catalog has been ex- 
tracted and compiled by S. Berta for the PEP collaboration 
(used in papers such as IWuvts et alj201 UlNordon et al.l2012l 
and others) and has been publicly released B T he GOODS-S 
optica l photometry is from the MUSIC catalog (ISantini et alJ 
2009). We exclude bands that have their central wavelengths 
between 2000-2300 A in order to avoid the bands from being 
dominated by an enhanced attenuati on feature around 2170 A 
(the so called 'UV bump ') that may dBuat et al.ll201 ll) or may 
not dCalzetfi et all 1 19941) exist at z > 1. To derive M\ m we 
use all remaining filters between rest 1200-2800 A and inter- 
polate the flux at 1600 A (converted to AB magnitudes) by fit- 
ting a power-law. The fitting minimizes a weighted^ 2 , where 
we weight the magnitude errors by 1 + |/tfii te r - 1600|/1600 A. 
This is done in order to give more weight to filters that ob- 
serve close to the desired wavelength. Errors are estimated 
by randomizing the input fluxes according to their respective 
errors, while assuming a normal distribution and interpolating 
again to obtain Mi6oo- The process is repeated a large number 
of times and the standard deviation on the derived M1600 is 
adopted as the error. To derive M2800 we select the filters that 
observe rest 1500-3500 A (again excluding the range of the 
'UV bump') and interpolate in the same way as we did with 
^1600- To ensure a reliable interpolation to the two wave- 
lengths we require at least one band that observes less than a 
full filter width from the desired rest frame wavelength and at 
least two bands in total. For this reason we limit our lowest 
redshift to z > 1, where 1600 A rest is covered by the U-band 
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Fig. 1 . — Several properties of the sample: left: The redshift distribution. 
Center. Stellar mass distribution. Right: Total infrared luminosity distribu- 
tion. 

and above. Finally, we select only galaxies for which the final 
error on /3 is lower than 0.2. 

Since the presence of an active galactic nuclear (AGN) can 
significantly affect the emission in UV, we remove galaxies 
suspected as AGN hosts from the sample. Sources which 
have an X-ray counterpart in the Cha ndra 2 Ms catalogs 
(Alexander et al J 120031: iLuo et al.l 12008) have conservatively 
been flagged as suspected AGN hosts and remo ved. None 
of the removed galaxies clearly comply with the Bau er et al] 
(2004) classification as star forming galaxies according to 
their X-ray emission. While we may be removing some ex- 
ceptionally rapid star forming galaxies, whose X-ray emis- 
sion is due to the st ar formation and not the AGN (e.g., 
lAlexander et al.ll2005l) . it will have a negligible effect on this 
study that focuses more on the bulk of the galaxy population 
around the main-sequence and a little above it. AGNs that are 
obscured enough to prevent detectable X-ray emission will 
also have their UV emission highly attenuated and the ob- 
served UV can be attributed to the star formation. 

In total, our sample includes 250 and 200 sources in 
GOODS-N and GOODS-S respectively. The properties of our 
sample in redshift, mass and Lir distributions are presented in 
figure Q] 

3. THE EFFECTIVE UV ATTENUATION VERSUS UV SPECTRAL 
SLOPE 

3.1. Terms and Definitions 

Star formation rates can be calculated both from the total 8- 
1000 jt/m infrared luminosity and from the UV luminosity at 
1600 A. We adopt the conversions of lKennicuttl(ll998 ). mod- 



ified downward by a factor of 1 .6 to match a lChabriei d2003l) 
initial mass function (IMF): 



SFR 



IR 



M yi- 



1.09 x 10 



_!0 MR, (8-1000 /jm) 



- = 0.875 x 10" 28 - 
Moyr- 1 



^v, 1600 A 



ergs 



-1 



Hz 



(1) 



(2) 



The conversion from Lir to SFR assumes that Lir represents 
the bolometric luminosity of the young stars. In practice, we 
detect UV emission, indicating that some radiation escapes 
and is not reprocessed into the infrared. Therefore, for the 
total SFR we use: SFR tot . dl = SFR lR + SFR VV . 

S FRuv uses the observed UV luminosity and hence is 
missing a significant fraction of the luminosity which was 
attenuated by the dust. The effective UV attenuation Airx 
should correct SFRuv to SFR tota i, i.e. log(5' FR totSL \) = 
log(SFR vv ) + 0.4Air X . We define Airx as: 



3 http://www.mpe .mpg .de/ir/Research/PEP/public_data_releases .php 
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Fig. 2. — The effective UV attenuation versus UV spectral slope for 1.0 < 
Z < 2.5 PACS detected galaxies. The color scale indicates the A log(SFR)MS 
distance of the galaxy from the main-sequence. L ines indicate various Afrx- 
B rela tions: solid black - Meurer et al. 1 1999), solid red - Cal zetti et al.l 
120001) . dashed green - best fit for galaxies with -0.15 < Alog(SFR) M s < 
0.15. The hatched area indicates expected values of/? for an unobscured pop- 
ulation. The lower panel shows the distribution in /? of the PACS detected 
sample (blue). The distribution of sources with similar selection criteria, but 
without a PACS detection is plotted in green. 

Let f> be the spectral slope, such that fx 00 AP. While spectra 
between two wavelengths do not always resemble a perfect 
power law, (3 can always be formally defined between two 
wavelengths as: 



toe(fM/fM) -0.4(Mi-M 2 ) 



-2 



(4) 



log(Ai/A 2 ) \og(AJA 2 ) 

where Mi and M 2 are the AB magnitudes at wavelengths A\ 
and Ai respectively. In this study we adopt A\ = 1600 A, A 2 = 
2800 A and obtain rest frame M\ and M 2 from the available 
photometry as explained in Section [2] 

3.2. AiRx-yS for The PACS Selected Sample 

In Figure |2] we plot the Airx-/? relation for the full 1 < z < 
2.5 PACS-detected sample. There is a significant scatter in 
the relation. However, the galaxies were color-coded accord- 
ing to their distance from the main-sequence in A log(SFR)MS 
and the color scale reveals a structure to this scatter. We define 
Alog(SFR) MS = log(SFR galaxy /SFR MS ), where SFR MS (M*,z) 
is the SFR on the main sequence for the given mass and red- 
shift. We use the Rodigh iero et al.l (1 20101) redshift-dependent 
main-sequence curves. These were derived from Herschel- 
PACS data in GOODS-N, partial to the data set used here. 
For convenience, we give the main sequence curves below: 
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Fig. 3. — The UV attenuation from IRX relation versus the best combination 
of /? and SSFR or A log(SFR)MS ■ Top: Best fit for the entire sample with no 
additional parameters. Center: When including SSFR as a parameter in the 
fit. Bottom: When including A log(SFR)MS as a parameter in the fit. The 
standard deviation in attenuation magnitudes from the 1 : 1 relation is indicated 
in each panel. 

To find the main sequence at a specific redshift we linearly 
interpolate between the above curves for the given mass. For 
1.75 < z < 2.5 we assume that the main sequence does not 
evolve beyond z = 1.75. 

It is evident that the higher the A log(SFR)MS of the galaxy, 
the higher it tends to lie on the Amx-/3 relation. As men- 
tioned in the introduction, local galaxies show a very simi- 
lar effect in the 'birth-rate' parameter, closely related to the 
SSFR. The main difference between SSFR dependency and 
A log(SFR)MS dependency is that constant SSFR is indepen- 
dent of the galaxy mass, while constant A log(SFR)MS is mass 
dependent due to the non-zero slope of the main sequence 
(e.q. [5}. To verify that indeed the slope of the main-sequence 
has an effect and that A log(SFR)MS is the more important pa- 
rameter for the scatter, we fit a plane to the data in Airx, /3 and 
either SSFR or A log(SFR)MS space. The result is plotted in 
Figure [3] As a first order correction, A log(SFR)MS allows a 
better minimization of the scatter than simply using the SSFR. 
Most of the effect is due to the slope and not so much due to 
the absolute scal e of the main-sequence, since the adopted 
main-sequence of Rodigh iero et alj (120101) evolves little from 
z = 1 to z = 2, but has a significant slope of ~-0.5 in SSFR 
versus mass. We warn that the relations in Figure[3]should not 
be used to derive Airx- Such attempts will result in a signif- 
icant fraction of large attenuation correction errors. We also 
expect systematic deviations, in particular for galaxies with 
j6 < -1.5 due to the dependency on A log(SFR)MS that may 
be non-linear over the full range of /?. 

Galaxies on the main-sequence form a near straight se- 



The FIR-UV relation up to z=2.5 



5 



quence in Airx-/? (Figure [2]). We fit a linear relation to the 
galaxies close to the main-sequence (-0.15 <Alog(SFR)MS< 
0.15) while minimizing the x 2 for the orthogonal projection 
of the Airx and ft errors and derive: 



A MS = 1.26/3 + 3.90 



(6) 



This slope is shallower than the relations derived by M99 
and COO (Figure [2]). Taken at face value and assuming that 
the obscuring material is arranged as a screen in front of the 
sources, an Airx-/? slope of o-trx = 1 .26 translates into a dust- 
extinction law with /T 1 18 wavelength dependency (e.q. lEll in 
the Appendix). Such a UV extinction law is almost iden- 
tical to the Small M agellanic Cloud (SMC) extinction law 
dGordon et al J 120031) . that has A(A) oc A~ 115 between 1600 
and 2800 A. 

It is clear that the way the sample is selected can signifi- 
cantly affect the derived Airx-/? relation. For example, a high 
FIR luminosity cut will tend to select galaxies above the main- 
sequence (in SSFR), resulting in a higher and possibly steeper 
Airx-/? relation. Mixing different redshifts can result in dif- 
ferent Airx-/? relation since galaxies with ULIRG-like IR lu- 
minosities have high A log(SFR)MS at z~0, but can be main 
sequence galaxies at z > 1 . Airx-/? samples that are selected 
to span a large range in ft will tend to select fundamentally 
different galaxies at the two extremes of the ft range. 

We also tested a subsample of BzK selected (Dad di et al.1 
2004) galaxies from our full sample. The picture for the BzKs 
does not change significantly from Figure [2] This is some- 
what expected because almost all z > 1 .4 galaxies in our sam- 
ple are also selected as BzKs, when the Ks imaging is deep 
enough. The M99, COO and Ams (ft) relations cross each other 
around ft ~ -1, which is also the peak of the number dis- 
tribution in ft (Figure [2] lower panel). The peak in the ft 
number distribution for sources undetected by PACS is at a 
slightly lower ft. This means that for the large majority of 
UV galaxies, all 3 relations will produce similar attenuation 
corrections. However, biases will develop depending on the 
assumed Airx-/? relation and the sample selection, in particu- 
lar for galaxies with much higher/lower ft than -1. 

In the next sections, we will study in more details the scat- 
ter in the Airx-/? relation, the parameters governing it, and 
the information it provides about the internal structure of the 
galaxies. 

4. MOLECULAR GAS AND THE SPECIFIC ATTENUATION 

Stars form in molecular clouds and a tight relation is estab- 
lishe d between the SFR and the mo lecular mass surface den- 
sity dSchmidllll9"59t lKennicutCl989l) . often referred to as the 
'Kennicutt-Schmidt' (KS) relation. Molecular clouds contain 
dust that is mixed with the gas and makes them optically thick 
to UV, obscuring the stars forming in it. The relation between 
the optical depth and the molecular gas content (with the dust 
embedded in it) strongly depends on the geometrical arrange- 
ment of the stars and dust. Geometries will be discussed in de- 
tail in Section [6] but for the moment we can consider the fol- 
lowing simple example: (a) two stars and two dust elements, 
each star with a dust element in front of it that produce a cer- 
tain attenuation, versus (b) two stars crowded behind a single 
dust element that obscures both to the same attenuation value 
as the former case. Therefore, (a) requires more gas and dust 
per obscured star than (b) to produce the same attenuation. 

Let us define a new quantity, the 'specific attenuation' (Sa) 
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Fig. 4. — The AirxtS relation for the CO sample . Colors and symb ols indi- 
cate the literatur e source: blue diamo nds are from[t)addi et al. 1 2010) , red tri- 
angles are from lTacconi et alj [ZMftk green squares are from[Magnelli et al. 
12012). The dotted line mark the Meurer et al. ( 1999) relation. The solid 
black line indicate the relation for the main-sequence galaxies from this work 
and dotted lines are spaced 1 magnitude difference from the main-sequence 
relation. 
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Sa is proportional to the attenuation per unit molecular gas 
available per UV emitting star. However, the total number 
of UV sources Nuv is n °t a convenient quantity to use. The 
calibration of the UV luminosity as a SFR indicator relies 
on the assumption that the UV emitting population reaches 
a stea dy state in less t han 10 8 years of continuous star forma- 
tion (iKennicuttll 19981) . Hence, the SFR in the above equation 
serves as a more practical proxy to the number of UV emit- 
ting stars (Nuv)- Tuv is a constant with time dimensions that 
makes Sa dimensionless. The exact value of Toy is of little 
consequence in this study as it will be blended by subsequent 
empirical fitting. Purely for scaling convenience, we will use 
Tuv = 10 9 yr. Thus, for the typical SFR/M mo , ~ 10~ 9 yr 1 
and Airx ~ 1, we expect Sa in order of unity. Sa tells us how 
efficiently the gas (and the dust embedded in it) is arranged in 
order to block the light from the stars. 

Sa strongly depends on the typical stars-dust geometrical 
arrangement (see Section|6]l and also varies with compactness. 
If we keep the number of young stars and amount of dust con- 
stant, and increase the local surface densities by compressing 
the dimensions of the star forming regions, then the column 
density will increase, resulting in a higher A/«x and a higher 
Sa- The physical meaning is that the stars are more crowded 
behind every gas (dust) element and each gas element blocks 
the lines of sight to more UV sources. Whether the stars are 
behind the gas or mixed with it will also have a strong effect 
on Sa- 

To derive Sa we require CO observation of galaxies at z > 
1 , that are also detected in the FIR. For the moment, we would 
like to avoid the extreme and rare cases of sub-millimeter 
galaxies (SMGs) and concentrate on more common galaxies 
near the main-sequence. A handful of such PACS-detected 
galaxies have CO measurements from the Plateau de Bure In- 
terferometer (PdBI) in GOODS-N and the Extended Groth 
Strip (EGS) fields. Table [2] specifies the details of our CO 
sample. The sources in GOODS-N were already included in 
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our m ain Airx-/8 sample and here we add the iTacconi et alj 
(2010) sources in the EGS field, also observ ed with Herschel- 
PACS as part of the PEP project dLutz et al.ll201 ll) . We correct 
the literature M mo i to use a uniform conversion factor orco = 
4.35 M Q K -1 krrT 1 s pc" 2 , such that M mo \ = aco ■ L co t_o, 
where Lco t-o ls m K km s pc 2 and M mo i is in M Q . aco al- 
ready includes a factor of 1.36 that accounts for the He mass. 
The optical photometry that is used in deriving j3 and 1600 A 
flux (same method as described in Section|2|i is from the pub- 
lic catalogs of the AEGIS0 team. 

In Figure|4]we plot the Airx-/? diagram for the CO sample. 
As can be seen, thi s sample is widely scattered around the 
iMeurer et~ai1 0999) relation and the A M s(fi) (e.q. [6jl relation 
discussed above. The sample has yfif > —1.5 for all the galaxies 
and avoids the cases of blue obscured galaxies, that suggest a 
bimodal distribution in /,),(A) (Section^. 

In section [3] and figure [2] we established that the loca- 
tion on the Airx-j6 diagra m is shifted upw a rds in A irx with 
increasing Alog(SFR)MS- ISaintonge et alj d201 lbl) found a 
relatively tight correlation between the gas depletion time- 
scale Tdep - M mo i/SFR and SSFR for local galaxies. This 
relation scales with redshift like the main-sequence, which 
means that more generally, a decreasing Tden correl ates with 
an increasing Alog(SFR) M s (Saintonge et al.ll20l2h . In fig- 
ure [5] top panel we plot SFR/M mo ; (t^L) as a function of the 

difference between the measured Airx (eq- [3]) and Ams(j6) 
(e.q. [6J. A clear correlation exists, which reflects the fact that 
A log(SFR)MS and Tdep are we H correlated. 

In the bottom panel of Figure|5]we plot S a as a function of 
Airx - Ams(j6)- A tightening of the relation relative to the top 
panel is evident. Note that although we are using Ams, only 
the slope in Ams is important since the zero order terms are 
adjusted when fitting a linear relation to Figure[5] The best-fit 
linear relation is: 



logOS*) = 0.20 ±0 .o4[Airx - 1.2603 - A))] + 0.46 ±0 
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(8) 



where fio = -2.2. The term A screen = 1.2606 — /?o) in the 
above equation presents 5 a as a function of the Airx excess 
over the attenuation of a uniform screen model for the given 
p. The rea sons for presenting S a in such a way are discussed 
in Section ROl 

While the correlation appears to be tight, the values in the 
x and y axes are correlated. In the bottom panel in particu- 
lar, Airx explicitly appears on both axes. To verify that the 
tighter relation in the bottom panel versus the top panel is in- 
deed significant, we perform a x 1 test that takes correlations 
into account (denoted^ 2 hereafter). A full description of the 
test is given in Appendix lAl The \\ values are indicated in 
each panel and indeed S a offers a significant improvement 
in tightening the relation. The x\ for the Sa relation is con- 
sistent with the expected from pure random noise of the flux 
measurements. Such a tight relation is non trivial. There are 
many possible arrangements of the stars and dust that will pro- 
duce very different Sa, while still having the same location in 
the AiRx-yS plane, yet this is not the case in our sample. 

The increase in Sa with increasing Airx at a given /? in- 
dicates that the shift of the galaxies on the Airx-/? plane is 
mostly due to geometry and not the unattenuated spectral 
slope /?(). There is no compelling reason for ySo to be tightly 
correlated with Sa- Variation in the extinction law due to dif- 
ferent dust grains between galaxies can in principle affect both 
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Fig. 5. — Top: The SFE (TjJ ) versus the attenuation distance from the 

AiRX-f> main-sequence relation. Colors indicate the Airx for the galaxies 
matching the color bar at the top. Symbols indicate the source for the gas 
masses. The reduced perpendicular-^ 2 is indicated in the lower right corner. 
Bottom: The same as above for the specific attenuation Sa- 

Sa (through A I Mmoi) and Airx at a given /3, to produce a cor- 
relation. However, Sa in our sample changes by factor ~6, 
which is more than what reasonable changes to the UV ex- 
tinction law can produce. We conclude that the scatter in the 
Atrx-jS relation is dominated by the geometrical arrangement 
of dust and stars as reflected by the values of the specific at- 
tenuation Sa- Other factors considered such as ySo and the UV 
extinction law seem to play a less significant role. 

5. ESTIMATING GAS MASSES WITH UV AND FIR PHOTOMETRY 

5.1. The S a Method 

The relation obtained in figure [5] bottom (e.q. [8]) combined 
with the definition of S a (e.q.|7]), can be used to derive molec- 
ular gas masses from integrated UV and FIR fluxes: 



lo 8 (^Fr) ~ 0-20(Air X - 1 .26/?) + 0.09 

(9) 

In figure [6] we plot the predicted molecular gas mass from 
e.q. [9] versus the measured values from CO observations. 
E.q. [9] gives a good prediction to the molecular gas mass 
for the galaxies in our z ~ 1 sample, on which the relation 
was calibrated (blue circles). The scatter is only 0.12 dex, 
with no bias. However, it is much more interesting to test 
this method for estimating molecular gas masses on different 
samples that were not used in its calibration. For this pur- 
pose we use a z^O subsample from the COLD GASS project 
(ISaintonge et al.ll201 lal 12012). local ultra luminous infrared 
galaxies (ULIRGs) and two z > 1 SMGs in GOODS-N. 

We select the COLD GASS galaxies that have FUV and 
NUV fluxes from GALEX, and 60 & 100 yum fluxes from 



4 http : //aegis . ucolick . org/astronomers . html 



5 http : //www .mpa- garching .mpg . de/COLD_GASS/ 
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TABLE 2 

THE Z ~ 1 SAMPLE OF NORMAL STAR FORMING GALAXIES WITH CO MEASUREMENTS AND HERSCHEL-PACS FLUXES. 



ID 


Source f 


z 


/? 


Airx 
mag 


10 12 L G 


sfrL 

M G yr 1 


M* , 

mo! 

10 10 M o 


BzK-4171 


D10 


1.465 


0.44 ±0.11 


5.6±0.2 


0.92±0.06 


101+7 


9.3±2.2 


BzK-21000 


D10 


1.523 


-0.60±0.10 


5.1±0.2 


2.10±0.08 


231 ±9 


9.8±2.4 


BzK- 16000 


D10 


1.522 


0.06±0.06 


3.8±0.2 


0.73±0.11 


82±12 


7.1 + 1.7 


BzK-17999 


D10 


1.414 


0.76 ±0.12 


6.5±0.2 


4.11±0.38 


450±41 


7.6+1.8 


BzK-12591 


D10 


1.6 


-0.13±0.05 


4.7±0.2 


2.41±0.10 


267+11 


14.5±3.5 


EGS 13004291 


T10 


1.2 


-0.10±0.06 


5.5±0.2 


5.32±0.14 


585+15 


28.0±1.1 


EGS 12007881 


T10 


1.17 


-0.22±0.05 


2.4±0.2 


0.77±0.15 


94±16 


8.3±0.4 


EGS 13017614 


T10 


1.18 


0.72±0.13 


4.1±0.3 


0.74±0.15 


82+16 


9.3±0.7 


EGS 13004661 


T10 


1.19 


-0.10±0.1 


3.2±0.4 


0.49±0.16 


57+17 


2.4±0.5 


EGS 13003805 


T10 


1.23 


0.60±0.15 


5.6±0.2 


1.93±0.21 


212±23 


18.0±1.2 


EGS 13011439 


T10 


1.1 


0.23±0.09 


3.6±0.3 


0.53±0.12 


60±13 


4.5±1.0 


GN21172 


M12 


1.249 


-0.50±0.12 


4.0±0.1 


0.42±0.04 


47±4 


1.62±0.44 


GN17022 


M12 


1.224 


-0.84±0.09 


2.7±0.1 


0.31 ±0.04 


37±4 


2.95±0.68 


GN1464 


M12 


1.084 


-1.15±0.17 


2.2±0.1 


0.36±0.05 


45±5 


3.63±0.91 


GN5532 


M12 


1.021 


-1.19±0.17 


3.9±0.1 


0.54±0.03 


60±3 


2.88±0.95 


GN26843 


M12 


1.016 


-0.92±0.17 


3.9±0.1 


0.98±0.04 


110+3 


5.75+1.38 


GN13168 


M12 


1.016 


-1.01+0.17 


2.3±0.1 


0.41±0.03 


51+3 


2.88±0.68 


GN177 


M12 


1.17 


0.63±0.11 


3.2±0.2 


0.26±0.04 


30±4 


5.13±0.72 



f DIP - Dad di et al. (20101 TIP - Tacconi et al. (20101 M12 - Magnelli et al. f20121 

* from Herschel-PACS data of the PEP project (TIP sources) and the combined PEP & GOODS- 
Herschel data (DIP and S12 sources). 

* SFR tot = SFR(FIR) + SFR(UV) uncorrcctcd 
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Fig. 6. — The M mo i(S a) as derived from the specific attenuation (S a) rela- 
tion, versus the M mo i(CO) measured from CO observations. The solid line 
is the one to one relation. Blue circles are for the z ~ 1 SFG sample, green 
circles are for the z ~ COLD GASS SFG. Red colors are used for extreme 
galaxies such as the z>l SMGs (star symbols), the z ~ COLD GASS 'ma- 
jor mergers' (squares) and z ~ ULIRGs (triangles). The standard deviations 
for the z ~ and z ~ 1 samples include only the SFGs. 

IRAS. This sample includes a few galaxies which are la- 
beled as 'mergers' that will be considered separately from 
the other star forming galaxies (SFGs). The local U LIRGs 
are selected from the sample of Solom on et al.l (fl9 97j) cross- 
matched with the sample of iHowell et al.l (120101) that pro- 
vide GALEX fluxes. All the ULIRGs are treated as 'major 
mergers'. Full details regarding the ULIRGs, the SFG and 
COLD GASS subsamples, as well as the derivation of IR lu- 
minosities and UV slopes are given in Appendix|B] A discus- 
sion regarding the conversion of IR luminosity to SFR in the 
COLD GASS sample is provided in AppendixICl 

We apply e.q. [9] on the COLD GASS subsample with no 
additional tuning. The result is plotted in figure[6] As we can 
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Fig. 7. — Histograms of the Sa values in the different samples. Sa values 
are calculated from its definition using e.q. [7] 

see, e.q. [9] which was derived for z-l galaxies, also pro- 
duces good molecular gas mass estimates for the z = SFG 
(green circles), which are an order of magnitude lower in M mo i 
and SFR than the z-l galaxies. The standard deviation is 
only 0.16 dex, slightly larger than the scatter at z — 1, though 
on average, e.q.|9]tends to underestimate COLD GASS M mo \ 
by 0. 1 dex. There are several possible causes to the 0. 1 dex 
offset. Changes to fio, to the UV extinction curve, or to the 
dust-to-gas ratio can potentially explain a small shift in the 
Sa~Ai RX - P relation between z-l and z — 0, that will 
affect e.q. [5] While in the z > 1 sample we avoided using 
filters that may suffer enhanced attenuation from a 2170 A 
'UV bump', we could not do so with GALEX. The NUV fil- 
ter is very broad and overlaps with a possible UV bump often 
observed in local galaxies. This will tend to lower which 
in turn leads to underestimation of M mo /. In order to explain 
the 0. 1 dex underestimation seen in Figure [6] (3 will need to 
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change by ~0.4, meaning NUV lowered by ~0.25 mag due to 
the 'UV bump'. 

'Major mergers' were not included in the z — 1 sample 
on which the method was calibrated. The derived M mo i for 
the z = 'mergers' (COLD GASS and ULIRGs, empty 
red symbols) are systematically over estimated, by an or- 
der of magnitude. Two SMGs in GOODS-N have PACS 
fluxes and good optical (rest-UV) photometry: HDF169 and 
HDF242 at redshif ts 1. 22 and 2.49 The C O fluxes are from 
iFrayer et aD (120081) and llvison et all (1201 ll) respectively. Be- 
ing PACS sources in GOODS-N and at the suitable redshifts, 
these sources are also included in our AiRx-yS main sam- 
ple (Section The two SMGs are also likely to be prod- 
ucts of major mergers which is also refl ected in their ago ~ 
1 Mq K" 1 km" 1 s pc~ 2 conversion factors dTacconi et al.l2008t 
Magd is et al.ll201 UlMagnelli et alJl2012l) . similar to the local 
ULIRGs. The two SMGs are plotted in Figure [6] as red stars. 
Similar to the local mergers, the SMGs M mo i are over esti- 
mated by the Sa method, though only by ~0.3 dex. 

Galaxies that are undergoing major mergers appear to fol- 
low a different Sa-Airx-/3 relation than the tight relation we 
find for normal SFGs. This likely reflects their compact, 
intense, central star forming region that dominates the total 
SFR. If we compare the Sa values (using arbitrary tuv = 
10 9 yr M Q ') we find a median of 6.4 and 3.4 M~ l for the z=l 
and z=0 SFGs, versus medians of 17, 40 and 134 M Q l for the 
SMGs, COLD GASS mergers and ULIRGs. Figure □ shows 
the Sa histograms of the various populations. There is little to 
no overlap in the spread of Sa between the SFGs and mergers. 
The two SMGs in our sample seem to be very moderate ex- 
amples and we expect other SMGs to reach higher Sa values, 
more similar to local ULIRGs. The large difference in abso- 
lute Sa in addition to the deviation from the Sa-Airx~/3 rela- 
tion is indicative of a change in the typical star dust geome- 
try, beyond a simple change in the compactness (local surface 
density) of the star forming regions. For example, a transition 
from scattered star forming regions with variable obscuration 
to a central star forming region with dust and stars mix can 
lead to a dramatic increase in Sa- We discuss this further in 
Section|6] 

E.q. [9] can therefore be used for molecular gas mass esti- 
mation on a wide selection of SFG galaxies. We must cau- 
tion that the method has been derived and tested on galax- 
ies that are mostly SFG with (3 > -1.5. Passive galaxies 
and/or very compact major mergers high above the main- 
sequence, may follow different Sa scalings, as is the case for 
the COLD GASS mergers and ULIRGs in our sample. We 
further discuss the possible application to major mergers in 
the next section. 

5.2. Comparison With The Kennicut-Schmidt Relation 

Another method commonly used for indirect estimation of 
gas mass is inverting the KS relation (inverse-KS). In order to 
compare the Sa method with the inverse-KS method, we use 
the three subsamples used above: The z~l CO sample (Ta- 
blegj and the z~0 FIR-bright COLD GASS subsample (Ta- 
ble [3j split into normal (labeled 'SFGs'f] and merger galax- 
ies. For the invers e-KS method, we adopt the relations of 
lGenzeletal.ld2010l) : 

logf -\ = 1.171ogf^|) - 3.48 (10) 

iMoyr^Kpc- 2 / 6 \M p C - 2 / 

6 Excluding COLD GASS ID 3962 due to lack of a reliable radius mea- 
surement 
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Fig. 8. — The deviation of the molecular mass, as derived by various meth- 
ods, from the molecular mass as measured from CO. Triangles are for assum- 
ing a constant gas depletion time scale, squares are for the inverse KS method 
and circles are for the Sa method in this paper. The standard deviation for 
each method in each of the three samples is indicated in the legend. 

log( ) = UV log fe^j- 2.44 (11) 

^Weyr-'Kpc- 2 / S \M pc- 2 j 

were e.q. [10] is used for the normal galaxies and e.q. [TT] for 
the major mergers. For the galaxies surface area we adopt the 
half light radii from the corresponding literature sources of the 
samples. The half SFRs are calculated from the S FR(IR) + 
SFR(\JY) as described in Section O 

In Figure [8] we plot the deviation of the estimated M mo \ 
from the value derived from CO measurements in dex. 
The standard deviation (dex) for each subsample with each 
method is displayed in the legend. The three subsamples are 
treated separately so possible systematic offsets between the 
samples will not translate into an increase in the standard de- 
viation. Biases can be adjusted by calibration, but the scat- 
ter is more inherent to the method. For the samples used 
in this study the S a method gives molecular gas mass esti- 
mates with an accuracy as good or better than the inverse-KS 
method, when avoiding extreme cases of major mergers. This 
is achieved using only integrated rest frame UV and FIR pho- 
tometry without the need to resolve the galaxies. 

For 'major merger' galaxies, the S a method has a large 
systematic bias in the derived molecular masses versus the 
inverse-KS method. The latter treats SFG and major mergers 
separately, using two different relations. However, for the lo- 
cal mergers, the scatter around the systematic bias of the S a 
method is not larger (slightly smaller in our sample) than the 
scatter in the inverse-KS relation, which suggests that it may 
be possible to derive a separate calibration for major mergers, 
that removes the bias. For example, in our sample, lowering 
the mergers M mo \ by 0.8 dex will place the S a mass estimates 
nearly on top of the inverse-KS estimates, object by object. 
0.8 dex is the difference between using e.q.[10]and e.q.fTTlfor 
these galaxies. Since here the 'major mergers' include only 3 
objects we cannot come to any solid conclusion regarding the 
application of the S a method to such galaxies. A much larger 
sample will be required in order to do so. 

6. DISCUSSION: DUST OBSCURATION AND GEOMETRY 

6. 1 . Effective Attenuation and The Attenuation Distribution 

In the previous sections we found that the scatter in the 
Airx-/3 plane correlate with the distance from the main se- 
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quence (Figure [2]i and also with the specific attenuation Sa 
(Figure|5]l. What is the cause of the scatter and why does Sa 
correlates so well with Airx - Ams? The physical parameters 
involved are the dust properties (extinction curve), the stellar 
population (unattenuated spectral slope) and the geometrical 
arrangement of the dust and stars. We must also keep in mind 
that here we do not observe a single star-forming region, but 
instead integrate the light from the entire galaxy. Hence, Airx 
is an effective UV attenuation that may differ from the attenu- 
ation of each individual source in the galaxy. To better under- 
stand this, we examine in this section the effects of integration 
on the AiRx-yS relation. Full derivations for the equations be- 
low are given in Appendix ITJ1 

When integrating the luminosity over an entire galaxy, we 
include many individual point-like sources and each may have 
a different attenuation value. Let us denote the normalized 
distribution of attenuations of the individual sources in the 
galaxy as /^(A), i.e. the fraction of sources suffering atten- 
uation between A and A + dA at wavelength A\. Let us now 
consider a general representation for fa . We will assume that 
most of the UV sources are scattered inside a typical range of 
attenuation values between A„„„ andA, m „+AA. For simplicity, 
we will assume that the distribution is fiat inside this range. 



fa{A) = 



1/AA 




else 



(12) 



Due to the integration over all sources, the resulting effective 
attenuation is not very sensitive to the fine details of fa lt with 
a possible exception around A — * 0. Most distributions that 
are unimodal can be approximated as fa above. 

We can now obtain the solution for effective attenuation at 
wavelength ,li: 



A Al =A min -2.5 log 



'log(e)(l-10- a4AA ) N 



0.4AA 



(13) 



By our definitions and when using the same rest UV wave- 
length in A\ and for deriving Airx, A,i, and Airx are equiva- 
lent. The observed spectral slope /3 is then given by: 



l 



logORu) 



-0.4(l -^ y 2 )A m ,„+log 



1 - 10 



-0.4AA 



1 - 10 



-0.4AAR, 



(14) 

where y is the power of the UV attenuation law (A{A) oc A y ) 
and R\ 2 — ^1 and A2 are the two wavelengths between 

which the UV slope is measured, 1600 A and 2800 A respec- 
tively in this study. With equations 1131 and [T4l above we can 
plot the expected Airx-/? relation for various fa(A) distribu- 
tions. 

6.2. Archetypal Dust and Star Geometries 

Let us now consider some possible, basic arrangements for 
the stars and dust mix. In figure [9] we sketch three archetypal 
stars and dust arrangements and one combination. These are 
by no means all the possible arrangements, but are sufficient 
for demonstrating the general principle. On the right hand 
side of the figure, we plot the approximation for the resulting 
fa(A). For simplicity, but still maintaining the essence of the 
matter, in all cases fa is approximated to be a flat distribution 
between A„„„ and A„„„ + AA. 

The two most commonly assumed geometries are the uni- 
form obscuring screen (case a) and the even-mix (case b) of 
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Fig. 9. — A sketch of three basic geometrical arrangements for the gas/dust 
and stars, a) The obscuring screen, b) Dust and stars mix, c) variable screens. 
a+b is a combination of a mix with an obscuring screen. The attenuation 
distribution for each case is plotted on the right hand side and parametrized 
by A„,i„ and AA 

stars and dust. The Airx-/? relation for both cases can be ob- 
tained from the general solution (e.q.[T3land[T4l. The special 
solutions for each case can be found in Appendix[E] 

In Figure [10] top panel we plot the expected Airx-/? rela- 
tions for uniform screen and even-mix models with various 
attenuation laws, while fixing /3o = -2.2. For the unifor m ob - 
scuring screen, the Airx-/? relation is linear with /3 (e.q. lEU . 
The slope is determined solely by y, while fio determines the 
intercept with the /? axis. Variation in the SSFR will change 
the ratio between the young and old stellar populations that 
tend to have lower and higher /?o respectively, resulting in a 
horizontal shift in the Airx-/? plot. Variations in fio can thus 
create scatter in the Airx-/? relation that is correlated with 
the SSFR and Alog(SFR)MS- However, even main-sequence 
galaxies at z > 1 have very high SSFRs, comparable with 
intense local starbursts and their UV light is still dominated 
by the very young population. To explain the large scatter 
in Airx by fio variations alone (under the uniform screen as- 
sumption) requires very low /3 values, lower than those seen in 
young stellar clusters and outside the typical p redicted range 
of -2.6 < (5 < -2.0 (e.g., IWilkins et all2012l) . The relation 
that we derived from the main sequence galaxies in section [3] 
implies /?o = -3.0, that is an unusually low value. However, 
the direct translation of the intercept in the linear Airx-/? rela- 
tion into /?o is relevant only if the geometry is that of a uniform 
obscuring screen, which is disfavored in this study. 

The even-mix geometry (Figure [10] dashed lines) also fails 
to explain the data. The even-mix distribution can only pro- 
duce very blue (i.e. low) /?. The limit for the change in the 
UV spectral slope is: 



lim /3-Po = -y 

AA— >co 



(15) 



With typical values of y 1, this distribution can only pro- 
duce a narrow range of low ft values, lower than in the ma- 
jority of our sample. It is clear therefore, that a broad fa (A) 
distribution that starts from near zero attenuation does not fit 
our sample. The typical galaxies in the sample must have 
fa{A) distributions that vanish when A — » and that effec- 
tively means at least some foreground screen component. In 
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P 

Fig. 10. — Top: Airx-/3 relations for a uniform slab model (solid lines) 
and even-mix (dashed lines). The colors are for various values of the UV 
extinction law slope y. The values along the dashed line indicate the A max 
value for the even-mix model. The data is plotted in gray in the background. 
Bottom: Airx-/3 relations for various values of a flat attenuation distribution 
starting at A„„„ with width AA. Lines of constant A„„„ are plotted in solid 
red. Lines of constant AA are plotted in dashed red. The relation for main- 
sequence galaxies is plotted in dash-dot green. The data is plotted in gray in 
the background. 



other words, A„„„ > and AA > are required. 

We denote the combination of even-mix with an additional 
foreground screen as case a+b. Such a distribution could 
arise from young stars embedded in molecular clouds start- 
ing at a minimal depth, or an additional obscuring component 
in the diffuse interstellar medium (ISM). Perhaps just as rele- 
vant when dealing with integrated luminosities, is the variable 
screens geometry as sketched in figure[9]c. In case c, each star 
forming region has a different screen and the region-to-region 
variations are what produces the AA width of fa(A), as op- 
posed to case a+b where AA is the depth inside a slab. As 
long as the resulting fa (A) in cases c and a+b are similar, they 
will produce a similar Airx-/? relation. Adopting y = —1.18 
from A MS (e.q.|6] >, B = -2.2 (between iMeurer et alll!999t 
ICalzetti et a l. 2000), and variable A„„„ and AA values, we plot 
the resulting Airx-/? relations in figure [TO] bottom panel. In 
this picture, z ~ 1 main-sequence galaxies have a typical at- 
tenuation distributions with full widths of AA ~ 12 magni- 
tudes, while galaxies with enhanced SFR have much broader 
distributions. The main sequence galaxies scatter along the 



linear relation mostly due to the A„„„ component that varies 
between them and controls the observed UV color. Being the 
most general case that combines both a screen and a broad 
distribution we can see in Figure[l0]that no A„„„ and AA com- 
bination can explain galaxies with B < —1.5 and Airx > 1. 
These galaxies either have a combination of lower Bo and a 
flatter attenuation law (higher y), or m ore li kely, a bimodal- 
like attenuation distributions (Appendix lE.3b . 

6.3. The Relation Between S a and Geometry 

On the Airx-/? plane there is no difference between a com- 
bination of stars-dust mix with an additional obscuring screen 
(case a+b) and the variable screens (case c), as long as the re- 
sulting attenuation distributions fa(A) are similar. There is 
however an important distinction between the Sa of the two 
cases: in case a+b, each layer of the slab obscure several stars, 
the observer-side layer obscuring the largest number. This al- 
lows a relatively small amount of gas to provide a high at- 
tenuation for most stars, as opposed to case c, where each 
star/region requires its own independent column of gas, re- 
sulting in a lower Sa (i.e., more gas per UV star to achieve 
the same integrated Airx). 

As we can see from e.q. [13] if the attenuation distribution is 
unimodal, then the deviation of the galaxy from the uniform 
screen (for which A screen = A„„„) Airx-/? relation behaves like 
~ log(AA), regardless if AA is an actual slab depth (case b) 
or a variation between regions (case c). The correlation of 
Sa with Airx - A screen (j8) is caused by the width of the fa (A) 
distribution, represented in our simplified model by AA, and 
regardless of A ,,„•„. An increase in AA move the galaxies away 
from A screen (/3) and corresponds to a more compact star form- 
ing regions or a stronger tendency towards a mix of stars and 
gas, thus Sa increases as well. Note that the values of A mm 
are significantly smaller than AA, and for most of our sample, 
A m ,„ < 2. This means that most of the gas mass is in the AA 
component and only a small fraction of the gas is in the fore- 
ground. This is likely the reason why Airx - A screen (j8), that is 
nearly proportional to log(AA), is such a good indicator to the 
gas mass in e.q. [9] 

The galaxies from which e.q. [8] was derived are all star 
forming disks, dotted with star forming regions and likely 
share a fairly similar basic stars-dust geometry. The same 
can be said about the COLD GASS SFGs. Moderate changes 
to the surface densities and/or geometry are enough to explain 
the limited range of Sa in these two SFGs samples (Figure|7|i. 
The local ULIRGs on the other hand are characterized by a 
central core of very dense star formation. The ULIRGs are 
also likely to have a fundamentally different obscuration ge- 
ometry, that results in their S a being higher by two orders of 
magnitude than that of the SFGs. In merger galaxies Sa in- 
crease both by the compactness of the star forming region as 
well as due to a radical change in the geometry type, thus they 
deviate from the relation found for SFGs (e.q. [8] and e.q. [9}. 
Mergers are also more likely to have a bimodal fa (A), which 
will throw off the relation between Sa and B. 

7. CONCLUSIONS 

In this work we studied the attenuation of UV (1600 A) 
light from star forming galaxies at redshifts 1 < z < 2.5 and 
its relation to the UV slope B. The effective attenuation Airx 
was derived from the FIR/UV ratio, also known as IRX. We 
use Herschel-PACS deep FIR imaging and optical photome- 
try from various instruments that cover rest-frame UV. Our 
sample shows a significant scatter in the Airx-/? relation, but 
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with a structure to the scatter, where galaxies with increas- 
ingly enhanced SFR relative to the main SFR-M* relation 
(main-sequence) show an enhanced Airx for their (3. 

When selecting only main-sequence galaxies, they tend to 
form a tighter sequence than the full population in the A IRX - 
P diagram, with a flatter slope t han commonly used rela - 
tions, such as lMeurer et all (1 1999b and lCalzetti et all (120001) . 
In terms of the UV extinction curve, our findings favor a 
steeper A(A) oc A 7 dependency than the above mentioned re- 
lations, with a power of y ~ -1.18, similar to that of the 
Small Magellanic Cloud. However, such a conversion from 
the Airx-/? slope to the extinction curve is geometry depen- 
dent. The above mentioned literature Airx-/? relations and 
the relation for main sequence galaxies found in this work 
cross each other close to the peak in the p number distribution 
of 1 < z < 2.5 galaxies. Hence, the attenuation corrections 
obtained from the different relations tend to produce similar 
global SFRs, but with different biases. 

Using a smaller sample of 1 < z < 1.5 galaxies with CO 
gas mass measurements, we studied the relation between UV 
attenuation and molecular gas masses. We find a very tight 
relation between the location in the Airx-/? diagram and Sa 
the 'specific attenuation', a quantity that represents the atten- 
uation contributed by the amount of molecular gas available 
per young star. Sa depends on both the geometrical stars- 
dust arrangement and the compactness of the individual star 



forming regions. 

We demonstrate that the tight empirical relation between Sa 
and the location in the Airx-/? diagram can be used to estimate 
molecular gas masses using UV and FIR photometry alone, 
without the need to resolve the galaxies. For normal SFG at 
< z < 1.5, the accuracy is between 0.12-0.16 dex when 
compared with masses derived from CO observations. Major 
mergers and SMGs seem to have a different stars and dust 
geometry and possibly follow a different S^-Airx-/? relation, 
that remains to be determined. 



PACS has been developed by a consortium of institutes 
led by MPE (Germany) and including UVIE (Austria); 
KUL, CSL, IMEC (Belgium); CEA, OAMP (France); MPIA 
(Germany); IFSI, OAP/OAT, OAA/CAISMI, LENS, SISSA 
(Italy); IAC (Spain). This development has been supported by 
the funding agencies BMVIT (Austria), ESA-PRODEX (Bel- 
gium), CEA/CNES (France), DLR (Germany), ASI (Italy), 
and CICYT/MCYT (Spain). 

This study makes use of data from AEGIS, a multi wave- 
length sky survey conducted with the Chandra, GALEX, Hub- 
ble, Keck, CFHT, MMT, Subaru, Palomar, Spitzer, VLA, and 
other telescopes and supported in part by the NSF, NASA, and 
the STFC. 
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APPENDIX 



A. THExl TEST 



This test uses a Monte-Carlo simulation in order to obtain the distribution of values which are derived from several observables 
with random errors. Both axes in Figure [5] are derived values with propagated errors and therefore must be treated equally. In 
addition, several observables are shared between the axes and therefore the propagated errors are correlated. By calculating x 1 
for the perpendicular deviation from a theoretical linear relation one can determine the 'goodness of fit' without bias towards 
any of the two variables (axes). By simulating a resampling of the observed data one can measure the distribution of the derived 
values from the simulation and do so simultaneously for both axes, deriving the distribution of the perpendicular distance from 
the tested theory. Thus, the correlations in the errors are also taken into account. 

We randomize all the input fluxes (UV, FIR, CO) according to their estimated errors while assuming a normal distribution. 
From the new set of fluxes we recalculate Tdep, Sa, Ams(/?) and Airx- We then calculate the perpendicular distance S±j from the 
original relations (in Figure[5]l for each galaxy i. The process is repeated a large number of times until we can derive a distribution 
for each S±j. The^ is then defined as: 



where 6±j is the perpendicular distance of the measured data point i from the mean relation and cr{6 ± j) is the standard deviation 
in 6±j as derived from the simulated distribution. For each simulated set of input fluxes we also re-fit the Sa (Airx - Ams) relation. 
The distributions in the values of the coefficients are fully consistent with the best fit and uncertainty values given in e.q.[8] 



We select the COLD GASS galaxies that have FUV and NUV fluxes from GALEX, and 60 & 100 //m fluxes from IRAS. 
Galaxy pairs, whose individual components are confused in some of the bands were excluded. f3 and the UV flux at 1600 A were 
calculated by fitting a power-law spectrum to the GALEX fluxes, while taking into account the (very broad) filters responses. 
The L/r(8 - 1000 /mi) luminosity was calculated by fitting CE01 templates to the two IRAS fluxes, while allowing both scale and 
template to vary. For a discussion on the conversion of Lir to SFR for this samp le see Appendix ICl Most of these galaxies tend 
to be SFG slightly above the z — main sequence (e.g., iBrinchmann et al.ll2004l) and we adopt aco - 4.35 M Kr 1 krrT 1 s pc~ 2 
conversion factor for these galaxies (similar to the z ~ 1 sample). A few of the selected galaxies are classified as 'mergers' and a 
different conversion factor {aco = 1.0 M Q Kr 1 krrT 1 s pc~ 2 ) was used in the derivation of their M mo /. The criteria for using the 
lower aco are IRAS flux ratio S(60)/S(100)>0.6 and Lir > 10 11 LO. Two of the three galaxies flagged as 'mergers' also show 
clear signs of an ongoing merger in their optical images. The essential information of the COLD GASS subsample is displayed 
in Tabled 

In order to include more ex treme local ga l axies, we build a sample of ULIRGs composed of objects shared between the samples 
of iSolomon etafl (fl997h and lHowell etall J20T0I). These ULIRGs were observed with GALEX and IRAS and we treat them in 
the same way as the COLD GASS sampleT TSolomon et al.l d 19971) found for these galaxies aco ~ 1.0 and we adopt this as a 
uniform value for all ULIRGs in the sample. See Table|4]for the list of objects and photometric fluxes. 

In Figure [TT] we plot the location of the COLD GASS and ULIRGs samples on the Airx-jS plane. The z~l sample is also 
plotted for c omp arison. Also plotted is the grid of Airx-jS relations for the various attenuation distributions, similar to Figure [TOl 



The IRAS selected COLD GASS sub-sample has a significantly lower SFR and SSFR (medians 7.7 M Q yr" 1 and 2 x 10~ 10 yr" 1 
respectively) than our z~l sample, although these values are still on the high side for local galaxies. The smaller young stellar 
population relative to evolved stars raises the question whether a sizable fraction the Ljr that we measure is due to radiation from 
old stars, reprocessed by the dust and emitted in the FIR. 

Several recent studies have attempted to address this issue and arrived at seemingl y contradicting resu lts, depending on the 
samples, wavelength range (FIR versus sub-millimeter) and methods. For example,[Smith et aQ (120 12h studied the SEDs of 
galaxies selected from the H- ATLAS survey that observe galaxies in He rschel-SPIRE sub-millimeter bands (250-500 /mi). They 
conclude that for galaxies with SSFR similar to o ur COLD-GASS sub sample as much as 50% of the 8-1000 /mi luminosity can 
be due to dust heated by old stellar populations. ICalzetti et al.l (120101) found an excess in 70 //m emission whe n comparing to 
the SFR as derived from a combination of Ha and 24 fxm emission, which they att ribute to heating by old stars. iLi et alJ (120101) 
studied the 70 /mi luminosity as a SFR indicator and by comparing to the results of Cal zetti et al.l (2010) conclude d that 25^4-0% 
of the total integrate d 70 /mi emission in their galaxies is not related to star formation. However, the samples of ICalzetti et al] 
(2010) and ILi etalJ d2010h were based on galaxies with typical SFR lower by an order of magnitude than the COLD GASS 
sub sample. If the luminosity due to star formation is sc aled up by an order of magnitud e, the fraction of 70 //m integrated 
emission due to old stars would go down to a few percents. lDommguez Sanchez et al.l d2012l) compared the SFR as derived from 
Ha and from Herschel-PACS 100 and 160 //m, for a sample of low redshift galaxies with very similar median SFR and SSFR as 
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TABLE 3 

The COLD GASS subsample I Saintonge et al. 201 1a) with GALEX and IRAS detections 



ID 


z 


NUV 


FUV 


P 


S60 


SI 00 




AlRX 


lo^ A/ 


a* 

CO 






AB 


AB 




Jy 


Jy 


10 10 L G 


mag 


M 




3819 


0.0453 


17.27 


17.89 


-0.54 


0.78 


1.44 


11.05 


2.66 


9.97 


4.35 


3962 


0.0427 


18.38 


19.09 


-0.34 


0.57 


1.39 


9.02 


3.7 


10.06 


4.35 


20286 


0.0346 


20.35 


20.71 


-1.15 


0.34 


0.99 


4.16 


4.96 


9.49 


4.35 


26221 


0.0317 


18.03 


18.81 


-0.18 


0.58 


1.8 


6.34 


3.69 


10.07 


4.35 


29842 


0.0341 


18.15 


19.09 


0.19 


0.27 


1.09 


4.37 


3.4 


9.68 


4.35 


24741 


0.0492 


17.64 


18.11 


-0.89 


0.25 


0.67 


5.26 


1.99 


9.62 


4.35 


6375 


0.0488 


17.46 


17.95 


-0.85 


0.38 


1.11 


9.57 


2.43 


9.78 


4.35 


6506 


0.0487 


18.1 


18.71 


-0.57 


0.32 


0.85 


6.57 


2.74 


9.99 


4.35 


40781 


0.048 


18.82 


19.22 


-1.06 


0.27 


0.82 


6.89 


3.32 


9.8 


4.35 


7493 


0.0264 


18.41 


19.12 


-0.34 


0.64 


1.53 


3.28 


3.7 


9.64 


4.35 


10019 


0.0308 


17.67 


18.15 


-0.87 


0.19 


0.7 


2.29 


2.15 


9.53 


4.35 


41969 


0.0351 


16.59 


16.9 


-1.27 


0.38 


0.93 


3.52 


1.31 


9.61 


4.35 


4216 


0.0438 


17.24 


17.64 


-1.06 


0.77 


1.96 


11.76 


2.58 


10.08 


4.35 


42013 


0.0368 


19.2 


20.25 


0.44 


0.75 


1.39 


6.94 


4.85 


9.96 


4.35 


4045 


0.0265 


17.52 


18.2 


-0.41 


1.01 


1.68 


4.07 


3.04 


9.71 


4.35 


19949 


0.0289 


18.67 


19.63 


0.24 


0.3 


1.61 


4.49 


4.31 


9.56 


4.35 


56632 


0.0281 


18.83 


19.81 


0.28 


0.23 


0.89 


2.4 


3.88 


9.5 


4.35 


14712 


0.0383 


17.14 


17.81 


-0.43 


0.32 


1.05 


5.45 


2.23 


9.82 


4.35 


18673 


0.0385 


18.69 


19.34 


-0.47 


0.38 


1.03 


4.79 


3.51 


9.59 


4.35 


50550 


0.035 


18.56 


19.11 


-0.71 


0.21 


0.73 


3.13 


3.06 


9.49 


4.35 


28365 


0.0321 


16.61 


16.95 


-1.2 


0.34 


1.15 


4.12 


1.62 


9.65 


4.35 


31775 


0.0412 


16.68 


17.23 


-0.71 


1.38 


2.71 


16.87 


2.68 


9.79 


4.35 


33039 


0.0434 


16.56 


16.89 


-1.22 


0.67 


1.34 


9.26 


1.74 


9.68 


4.35 


12533 


0.0331 


16.64 


17.17 


-0.75 


0.35 


1.16 


4.44 


1.79 


9.46 


4.35 


30332 


0.0425 


16.7 


17.12 


-1.01 


0.6 


1.27 


8.3 


1.86 


9.95 


4.35 


13005 


0.0487 


17.4 


17.96 


-0.68 


0.52 


0.97 


8.63 


2.34 


9.84 


4.35 


25327 


0.0314 


16.73 


17.27 


-0.73 


3.35 


5.07 


17.8 


3.34 


9.73 


1.0 


13666 


0.0405 


16.72 


17.19 


-0.89 


1.2 


1.84 


10.79 


2.25 


8.99 


1.0 


44942 


0.0329 


16.42 


16.88 


-0.92 


0.89 


1.86 


7.2 


2.0 


9.74 


4.35 


1115 


0.0259 


15.97 


16.4 


-0.99 


0.62 


1.6 


3.69 


1.47 


9.62 


4.35 


24973 


0.0285 


16.16 


16.99 


-0.06 


2.89 


5.59 


16.37 


3.16 


9.97 


4.35 


1221 


0.0345 


17.16 


17.83 


-0.43 


0.97 


1.63 


6.72 


2.66 


9.73 


4.35 


41904 


0.0392 


18.42 


19.12 


-0.36 


4.26 


3.49 


29.51 


5.18 


8.95 


1.0 



' H2+He mass. 

* The conversion factor from Leo l [K Km s _1 pc 2 ] to M„ M ; [M ], in M Q K~' kirT 1 s pc~ 2 

units. A value of 4.35 is adopted for the normal SFG and 1.0 for galaxies classified as 'mergers', 
aco includes a factor of 1.36 that accounts for the He mass in the H2 phase. 



TABLE 4 

The combined local ULIRGs sample of Solomon et al. 1 1997) and Howel l~et alJ 120101) . 



ID 


z 


NUV 


FUV 


p 


S60 


S100 




Airx 








AB 


AB 




Jy 


Jy 


10'°L G 


mag 


M G 


UGC 05101 


0.0394 


0.0 


0.0 


-1.05 


11.54 


20.23 


1.08 


6.78 


9.67 


IRAS F10565+2448 


0.0432 


0.0 


0.0 


0.11 


12.12 


15.13 


1.19 


7.49 


9.75 


UGC 08387 


0.0234 


0.0 


0.0 


-1.69 


15.44 


25.18 


0.48 


5.42 


9.58 


Arp220 


0.0182 


0.0 


0.0 


-0.93 


103.8 


112.4 


1.64 


8.05 


9.85 


NGC 6240 


0.0244 


0.0 


0.0 


-2.0 


22.68 


27.78 


0.69 


5.28 


9.9 


IRAS F23365+3604 


0.0646 


0.0 


0.0 


1.66 


7.09 


8.36 


1.53 


6.35 


9.93 



H2+He mass. Assumes aco = 10 



our COLD GASS subsample. They find an excellent agreement between the two SFR indicators, in particular for spiral galaxies 
like our COLD GASS galaxies. Since old stellar populations produce very little Ha, their result indicate that in their sample, very 
little of the Ljr as de rived from 100 and 160 /jm (very similar to the method used in this paper) is due to old stellar populations. 
Boqui en et alJ {2010) studied in detail the 100 and 160 /im emission in M33 and conclude that the 100 /im emission is a good 
tracer of SFR, while 160 fim emission may be affected by other processes. The origin of the FIR emission from region of low 
SFR density is unclear and the interpretation is strongly affected by uncertainties and modeling. 

Correcting the derived Ljr to remove emission which is unrelated to star formation is therefore not trivial. Studies based on 
FIR bands (60-100 //m) tend to find smalle r contamination from old stellar population than studies based on sub-millimeter 
(250-500 fim) bands. As Smi th et all (|2012) point out, sub-millimeter selected galaxies tend to be colder and dustier than FIR 
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O Z-1SFG 
• CG SFG 
■ CG Mergers 

A z-0 ULIRGs 

const. A min 

- - ■ const. A.4 



Fig. 11. — -Airx-/3 relation for the local samples: COLD GASS (CG) SFGs in blue circles, COLD GASS 'mergers' in red squares, local ULIRGs in mage nta 
triangles. The z~ 1 CO sample is plotted in empty green diamonds for reference. The gray grid indicate lines of constant A m j„ and AA, identical to Figure 10 



selected galaxies, and their L/« is more likely to be affected by radiation from old stars. Diffuse dust heated by the old stellar 
population is expected to be much colder than dust which is heated by the young stars and will dominate the emission only at long 
wavelengths, in the sub-millimeter bands. L/« as derived from SED fitting to the shorter FIR wavelengths, that are dominated by 
star formation related emission, does not include a component of cold diffuse dust. The fitted SED may under-predi ct the sub- 
millim eter luminosity, but still provide the correct L m to be converted to S FR using the bolometric consideration of lKennicuttl 
( 1998). Since for the COLD GASS subsample we fit lCharv & Elbazl d2001l) SEDs to the IRAS 60 and 100 fim bands and since 
the SFR and SSFR of these galaxies is relatively high for local galaxies, we do not apply any correction to the derived Ljr before 
converting it to SFR. 

D. DERIVATION OF AN ANALYTICAL Amx~P RELATION 

Let us define A A the effective attenuation at wavelength A as: 

L A = L A , sourc ,xlO-° AA < (Dl) 

A A therefore represents the magnitude by which the observed luminosity L A is attenuated with respect to the integrated luminosity 
at the source L AtSOllrce . The luminosities are integrated over the entire galaxy and A A includes all effects of absorption and 
scattering. Generally, A A will be wavelen gth d ependent and A A] + A^ 2 . As a result, observed flux ratios (colors) will differ from 
those ratios at the source. Combining e.q, IDll andl4l we can write: 

I -A, = -0.4 ^'"^ (D2) 

l0g(/li//i2) 

where j3 and fio are the observed and unattenuated UV spectral slopes, respectively. 

When integrating the luminosity over an entire galaxy, we include many individual point-like sources and each may have a 
different attenuation value. Let us denote the normalized distribution of attenuations of the individual sources in the galaxy as 
fa (A), i.e. the fraction of sources suffering attenuation between A and A + dA at wavelength A\ . We then apply e.q. lDll on each 
source in the galaxy according to the fa (A) distribution and integrate to derive the total observed (attenuated) luminosity. The 
effective attenuation that converts between the total emitted (unattenuated) luminosity and the total observed luminosity is: 

00 

A M = -2.5 log J fa(A') x 10-°- 4A 'dA' 

°oo (D3) 
A Al = -2.5 log / fa(A') x lO-O-^'/^dA' 

o 

Here we obtain A Al through the attenuation distribution in A\ with the use of the extinction law (represented by y) to convert 
from A' A2 to A', under the integral. Extinction laws in the UV-to-optical wavelengths can be approximated as a power law with 

power y, excluding a 'UV bump' around 2170 A in some examples. For a point source and any two wavelengths, y can always 
be defined through the relation: A Al jA Al - (A\ /A2) 7 , even if the extinction curve is not a true power-law. 

Let us now consider a general representation for fa . We will assume that most of the UV sources are scattered inside a typical 
range of attenuation values between A,„,„ and A„„„ + AA. For simplicity, we will assume that the distribution is flat inside this 
range. 

/il(A) - else (D4) 



The FIR-UV relation up to z=2.5 
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Due to the integration in e.q. lD3l the effective attenuation is not very sensitive to the fine details of fa , with a possible exception 
around A — > 0. Most distribution s tha t are not strictly multi-modal can be approximated as fa above. 
We can now insert fa into e.q. lD3l and obtain the solution for effective attenuation at wavelength A \ : 



■2.5 log 



fl og( e )(l-10- 
0.4AA 



(D5) 



The solution for A^ can be obtained simply by multiplying A,,,,-,, and AA in the right hand s ide o f the above equation by (Ai / A2) 1 ■ 
The observed spectral slope [5 is then obtained by inserting the above A^ and A^ into e.q, ID2l 



log(fii, 2 ) 



/ 1 - \Q-° AAA \ 

-0.4 (l - R~\)A min + log 

V 1,27 mm \i _ iQ-OAAAR^lj 



(D6) 



where R12 = A1/A2. By our definitions and when usin g th e s ame rest UV wavelength in A x and for deriving A IRX , the two 
attenuations A^ and Airx are equivalent. With equations |D5 1 and |D6l above we can plot the expected Airx-/? relation for various 
fa distributions, as we did in Figure [TOl 

E. A IRX -/3 RELATIONS FOR ARCHETYPAL CASES 
The Uniform Obscuring Screen 

The simplest case to consider and the one which is most commonly assumed, is a uniform obscuring screen (Figure |9j case 
a). In this case, all sources are subjected to the same attenuation magnitude and the distribution has the form of a Dirac delta 
function fa (A) = S(A^). We set A m ,„ = A^,, AA — > and the log term in e.q.[13]and e.q.[14]goes to zero. The Ajrx~/3 relation 
for wavelength Ai becomes: 

A~=- f yyff p-A) (ED 

1 - (AI/A2) 7 

Notice that the selection of the two wavelengths that are used in the definitions does have a small impact on this relation. The 
Airx-/? relation for this case is a straight line. The slope is determined solely by 7, and /3o determines the intercept with the /3 axis. 
We stress that the association of y with the linear slope and /?o with the intercept is specific to the obscuring screen geometry and 
may not be so in other cases. 

The Even-Mix 

Another case which is often considered is one in which the stars and dust are occupying the same slab volume with a constant 
densities ratio. It is represented by sketch b in Figure|9] In the special case of uniform densities, the attenuation distribution is: 

, f 1/AA 0<A<AA 

^^={o else (E2) 

When integrating the light over the entire galaxy, such a broad distribution makes some physical sense. Stars are forming inside 
molecular clouds at various depths inside the clouds, not only at the core. Older stars and stars that form near the edge of the 
cloud can clear away much of the obscuring material in their immediate vicinity and lower the attenuation they suffer and as a 
result, there is some width to the attenuation distribution. It is also relevant for galaxies that are observed nearly edge-on if the 
star forming regions are scattered in the disk and the inter stellar medium contributes a significant amount of attenuation. The 
expected result in these cases is a distribution of optical depths, starting from near zero and up to the maximum depths of the 
clouds, approximated as a step function in e.q, IE2l 
= 0ine.q.[H] 



-2.5 log 



'log(e)(l-10 



-0.4AA\\ 



0.4AA 



(E3) 



A similar solution can be obtained for A,\ 2 by replacing the AA terms with AA(A\IA2Y y m e.g. IE3I For AA >> 1 the effective 
attenuation increase like A^, ~ log(AA). The even-mix distribution can only produce very blue (i.e. low) ft. The /? value for such 
a distribution is given in e.q.[l4] For the limit A„„ ;l = and AA — > 00, we obtain: 

lim p - ft = -T (E4) 

AA— >oo 

The Bimodal Distribution 

Another simple case is a bimodal distribution fa = a6(A\) + (1 - a)6{A2). This is relevant for a case in which some of the 
stars are forming in large dense clouds, for example at the galactic center, and others in normal HII regions of low obscuration, 
scattered in the disk. 

A bimodal distribution has three free parameters, but since we have only two flux measurements there are an infinite number 
of solutions. For simplicity, we will choose the most trivial solution that has only two free parameters: 

fa = a5(A screen ) + (1 - a)6(A -> 00) (E5) 
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where a is the fraction of the population suffering the lower attenuation. This distribution has one component at an 'infinite' 
optical depth that does not contribute to the observed UV fluxes and the other at the correct optical depth to produce the observed 
jS. Here, 'infinite' means high enough so that the observed light is dominated by the component at lower depth, even if the 
obscured one is much more luminous at the source. 

While such a distribution cannot be represented in terms of A min and AA, it is easy to obtain the solutions for A Al and A Al from 
e.q.|D3] We get: 

A Al = A Ucieen (J3) - 2.5 log(a) (E6) 

where A scleen (yS) is the solution for the uniform screen geometry (e.q. lEU . Since a < 1, we get A A] > A scleen (yS) and if we attempt 
to derive the attenuation while assuming a uniform slab geometry we will underestimate the attenuation. On an A^x-ft plot, an 
increasing obscured fraction (decreasing a) will appear as a vertical shift in the attenuation for a given /? value. This distribution 
can produce any [} that the obscuring screen can produce and any increase in A above the A^x-P relation for the obscuring screen. 
Even very blue galaxies with [3 « fio can have significant attenuation and such a distribution can explain the galaxies with blue 
UV colors and relatively high attenuation (J3 < —1.5, Airx > 3). 



